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Radial slave boson representations have the particular advantage that the expectation values of their re- 
spective fields are finite even without the formal introduction of spurious Bose condensates for each of the 
bosonic fields. The expectation values of the radial (real) fields are in fact to be interpreted as the density 
of empty or singly occupied sites. Whereas the radial representation of the Barnes slave bosons has been 
investigated before, a setup for the functional integral of radial bosonic fields in the more physical Kotliar- 
Ruckenstein representation has not been accomplished to date. We implement a path integral procedure 
with suitable renormalization factors for a strongly correlated two-site model which allows to control the 
formal steps in the intricate evaluation, as the results for the partition function and the expectation values 
are known from exact diagonalization for such a minimal single impurity Anderson model. The partition 
function is shown to be a trace over a product of matrices local in time and therefore can be calculated ana- 
lytically. Eventually, we establish the scheme for the evaluation of correlation functions and thermodynamic 
properties. 

Copyright line will be provided by the publisher 



1 Introduction 

Application-oriented, fundamental properties of correlated electronic systems have become numerous in 
recent years, in particular for transition metal oxides: high-T c superconductivity (see, e.g., HED, colossal 
magnetoresistance (see, e.g., Q), transparent conducting oxides (see, e.g., B), high capacitance het- 
erostructures |]5] or large thermopower (see, e.g., [0), to quote a few. In addition, they also entail fasci- 
nating phenomena such as superconductivity at the interface of two insulators (7), peculiar magnetism in 
low dimensional systems [8], high temperature ferromagnetism in vanadate superlattices (9), all of them 
providing a strong challenge to investigate these systems from the theory side. Yet, it is fair to say that 
current theoretical approaches meet with severe difficulties when studying the models which describe these 
systems. Indeed, the tool, which is best mastered (perturbation theory), badly fails when the Coulomb in- 
teraction, both local and non-local, is sufficiently strong, thereby calling for alternative approaches. The 
introduction of generalized coherent states ifTUI in a constrained Fock space allows for the setup of a well- 
defined functional integral for strong Coulomb interaction ifTTIl but these path integrals are difficult to 
handle on account of the non-trivial Berry phases. An appealing scheme, which is rather based on canoni- 
cal fermionic and bosonic fields, the slave boson representation of local states, was introduced by Barnes to 
tackle the single impurity Anderson model (SIAM) 1121 . Later the slave boson approach has been extended 
to a whole series of models 01344211 (for a review see El ) . 



Corresponding author E-mail: Raymond.Fresard@ensicaen.fr, Phone: +33-(0)231 45 26 09, Fax: +33-(0)231 95 16 00 



Copyright line will be provided by the publisher 



2 



R. Fresard and T. Kopp: Exact results 



Notably, Kroha et al. 11231 developed a scheme which guarantees local gauge invariance in a conserving 
approximation and allows for Fermi liquid as well as non-Fermi liquid behavior for the investigated multi- 
channel Anderson impurity problem. In a different approach towards local gauge invariance, Read and 
Newns demonstrated that the phase of the slave boson field may be gauged away, at the price of introducing 
a time dependent constraint 11241 . By doing so the slave boson is reduced to a radial slave boson field. Such 
fields do not possess dynamics on their own, but only through their coupling to the fermions. The radial 
representation has not been intensively investigated, and the corresponding literature mostly focuses on the 
Barnes representation of the SIAM 112511261 . 

Besides, there seem to be (infinitely) many slave boson representations to a particular model. This has 
been shown to be of advantage by Kotliar and Ruckenstein when investigating the Hubbard model lfT4l . as 
they devised their representation in such a way that the Gutzwiller approximation is recovered in a para- 
magnetic saddle-point approximation. This approach is exact in the large degeneracy limit Q17II27I and 
obeys a variational principle in the limit of large spatial dimensions where the Gutzwiller approximation 
and the Gutzwiller wave function are identical 11281 . Extensions 1 1 6-18 20] have been applied successfully 
to a number of problems H291I30I . all of them making use of a specific representation of the kinetic energy 
operator. As the suggested operator takes an involved form, which is only defined through a corresponding 
series expansion, it could not yet be proven to comply to the normal ordering procedure, and therefore 
deserves further study. This is the aim of this paper. The correctness of various representations of the 
kinetic energy operator and interaction terms will be judged by the exact evaluation of the resulting func- 
tional integral representation of the partition function. We then illustrate the calculation of the expectation 
value of the radial slave boson fields. Being intrinsically phase stiff these fields do not average to zero, in 
contrast to canonical bosonic fields. Therefore their finite expectation value is by no means related to a 
Bose condensate. 

2 Gauge symmetry and radial slave boson representation 

We investigate an exactly soluble model, the simplified single impurity Anderson model enriched by a 
non-local Coulomb interaction term. The hamiltonian reads: 

U = ^2 (ec4 c CT + e & a i a a + v ( c t a a + h- c -)) + Ua \ a t a l a l + In ^ n c (1) 

rr 

which is defined on two sites. The one-body terms involve the operators cj. (c a ) and a\ (a a ). They repre- 
sent the creation (annihilation) of "band" and "impurity" electrons, with spin projection a. The energies 
e c and e a are the band and impurity energy levels, respectively, while the hybridization energy is V. Re- 
garding the two-body terms, / represents the non-local Coulomb interaction and U is the on-site repulsion. 
It is the largest energy scale in the model and it will be set to infinity as in the standard infinite-?/ SIAM. 
Standard perturbation theory in the interaction may not be considered for this problem. While this model 
may be solved by means of exact diagonalization ET1 . we here investigate it by means of a Kotliar and 
Ruckenstein slave boson representation. There, the local physical electron operator a a for the impurity 
site is rewritten as a product of the auxiliary fermionic fields f a and auxiliary bosonic fields e, p a and d, 
representing the empty, singly occupied and the doubly occupied site, respectively. The partition function 
reads for U — > oo: 



x e - (2) 
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with the Lagrangian 

£(r) = C f (r) + C b (r) + C v {r) + d(r) (3) 
consisting of 

£f(r) = £/j(T)(d r -,i + z:A( 2 ))./- CT (T) 

C b {r) = et(T)(9 r + zA( 1 ))e(r)+^pt (T) ^ T+iA (i) - U^) p» - *A« 

a 

C v {t) = y^( e t(rK(T)/ CT (r)4(r)+ e (r)pt( r)/ t (r)Ca(T) ) 

a 

d(T) = /(l-et(r)e(T))^4(r) C(T (T). (4) 

cr 

The first two contributions entail the dynamics of the auxiliary fields, together with the constraints, spe- 
cific to the Kotliar and Ruckenstein setup. The hybridization and the non-local Coulomb interaction are 
represented by the last two terms. Here, the slave bosons e {p a ) refer to empty (singly occupied) sites, 
respectively. They are subject to the constraints 



fJa-PlPa = <7=t4 • (5) 

They are enforced by the Lagrange multipliers A^, resp. A^ 2) . The La grangian, Eq. (01, possesses a 
U(l) x U(l) x U (1) gauge symmetry group. Indeed, performing the gauge transformation of the fermionic 
fields 



/» — ► e-^Ur) (6) 
and expressing the bosonic fields in amplitude and phase variables as 



e(r) = ^Rjr)e W[T) 
p a (T) = ^RAr)e^^ +9 ^ (7) 

allows to gauge away the phases of the three slave boson fields provided one introduces the three time 
dependent constraints: 

a(r) ee i (A (1) + <9 t (9(t)) 

Mr) ee -drxArj) ■ (8) 

We here introduced the radial slave boson fields in the continuum limit following, e. g., Refs. Il24l . When 
dealing with discretized time steps special care has to be taken. Extending the procedure introduced in 
Ref. [25] for Barnes slave bosons to the Kotliar and Ruckenstein representation one obtains the partition 
function as: 
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Z = lim lim lim (f[ F ^ F **jt* F ^ f° di?e „ f° ^ /°° dR J 



x e~ bb Z f 

where i? e „ {R a , n ) corresponds to the amplitude of the above complex e„ (p a ,n) bosonic variable, 5 
0/N, and' 



(9) 



N 



Sb 



(-R e ,n + -Rf ,n + ^4-,n — 1) — PanRa,' 

a 

l[D[c a ,ci]D[f a Jl]e- s f. 

a 

The fermionic contribution to the action is bi-linear in the fermionic fields. It is given by 



(10) 



N 



71=1 ET 

+<4 jff [c^ff - Cn_i )(T [l - (5 (e c - ix + 1(1 - y 2 a 



(ID 



Hence the fermions not only experience fluctuating magnetic fields, but the hybridization is effectively 
fluctuating as well. As several forms for the operators y an and z an may be considered we postpone their 
actual definition. While there is little intuition linked with the evaluation of Eqs. (l9l-(TTTTi the partition 
function may be cast in the particularly suggestive form 



Z = lim lim lim VZt 

Sl->oo N^oo n— >0+ 

namely as a projection of the auxiliary fermion partition function. Here we introduced 



(12) 



(13) 



with 



V n = 



f" ^ pug). r dR ^ r dRui r dRu , 

J-n 2lT J-Q 27r J-Q 27T J-n J-v J-v 



x e 



»<5(a„(fle,n+«t,n+fl4.,»-l)-E„ PonRa.n) 



(14) 



which is defined on one time step only. Having written the projection operator in the simple form Eq. (fT3t 
helps to evaluate the partition function Eq. (fT2l . Nevertheless this may not be straightforwardly achieved 
as the auxiliary fermion partition function consists of all up-spin and down-spin world lines describing 



Copyright line will be provided by the publisher 



andp header will be provided by the publisher 



5 



fermions in fluctuating fields. Yet, in Ref. Il26l . it has been shown that all these world lines may be re- 
summed when expressing Zf in the form: 



JV \ / N \ / N 



Zf = Tr J] K t „ • Tr JJ JC in = Tr JJ /C t „ ® IC in (15) 



n— 1 / \ n— 1 / \n— 1 



where the time steps are apparently not mixed. The re-writing of the fermionic determinant Eq. (fT5t from 
the action in Eq. (fTTT > may be gained by direct calculation and yields a convenient form to perform the 
exact projection onto the physical Hilbert space. The matrices JC an read: 



fcan 



/ 1 \ 

A an 5Vz an 
5Vz an L 

an 



(16) 



We introduced the short-hand notations: 



L an = e-^-^M = L a e- iS ?™. (17) 
At this stage we are ready to test various forms that may be taken by the fields z an , z an , and y a 



first concentrate on z an and z an . In the standard Barnes representation they are both given by R e , n , even 
though, in Cartesian gauge, they should be replaced correctly by (translated by e* ) and e (translated by 
e n _ 1 ), respectively (see Ref. Il25l ). In the Kotliar and Ruckenstein representation the operator form for z a 
is z a = e^p a + p_ a d (d — s- in the limit U —> oo). In accordance with the above result for the Barnes 
representation one would naturally suggest: 



Zan = \J Re,nRa,n+l (18) 

This forms respects the assignment that when an electron hops onto the impurity its occupancy changes 
from empty (y/R e , n ) to singly occupied (yj Ra, n +i), while when the electron leaves the impurity its occu- 
pancy changes from singly occupied (y/ R a , n -i) to empty (y/R e . n )- Of course, the assumption Eq. ( TT8l 
needs to be verified explicitly, for instance by computing the partition function Eqs. (fT2HT5l l. However an 
exact evaluation turns cumbersome since several time steps are involved in the matrices JC an and a large 
number of world lines need to be summed up. Yet, second order perturbation theory in V may be carried 
out, with the result that the correct answer is in fact obtained with the assignments Eq. ( TT8b up to terms that 
vanish in the thermodynamic limit. Though encouraging, this result does not prove the correctness of the 
representation Eq. (fT8l and, at least, higher orders in V should be considered. But, as infinite order seems 
to be beyond reach, we here follow another route. 

Let us note that, in their original paper, Kotliar and Ruckenstein used a different form for the z-operator. 
In operator form, for finite U, it reads: 



e f (1 -pXp„ - dU ) 2 (l-eh -pl a p_ a j 

ptril-plPr-rfdy* (l-e+e 2 d (19) 
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In the limit U —> oo, d — s- 0. Introducing radial slave boson fields, Eq. ([T9l — jointly with the corresponding 
expression for z\ — may be transformed into: 



z an — z an — \ ~j (20) 

y i — n a ,n 

Regarding the hole density that is represented by y^ n , a natural translation is: 



Van = y/Re,n (21) 

With this form the time steps are neither mixed in the projection operator Eq. (foi l nor in the fermion 
partition function Eq. 1151 . Hence one may write the partition function as 



Z = Ti(k N ) (22) 

with k = V n • ^nf ® fcn± where k is n-independent, respecting the translational invariance in time of the 
model. Therefore we need to project the entries of this matrix. Explicitly, the action of this projection on 
the various contributions entering Eq. ([TBI are found to be: 

' Zan 1 7~*n ' ^~<rn^—an -^c 7~*n ' ^an^—an (23) 

' Van 1 ^Vi ' ^trn^—an L/ c L a G l~*n ' L an L — an 

"P. A - T T> . A A T - T 2 T P - 2oJ T> . J A - T J p~ SI 

where L c , L a , and L an are defined through the relations in Eq. (TTTt . We obtain 



k = ( 1 ) ® ( xv I ) ® ( % I ) ® ( L ^~ SI ) © ( L ^~ 51 ) 

L^e" 2 " ) (24) 



where the direct sum (©) relates to the block-diagonal form of k that follows from labeling its rows accord- 
ing to: (l t , 1 4 ) -> 1, (l tl 2 4 ) -> 2, (l t , 3 4 ) -> 3, (l tl 4 4 ) -> 6, (2 t , 1;) -> 4, (2 t) 2j.) -> 8, (2 t , 3;) -> 
9, (2 t ,4 4 ) ->■ 11, (3 t ,l ; ) ->■ 5, (3 t ,2j 10, (4 t , l x ) 7, (4 t ,2 ; ) 12, with z T labeling the 
rows of /C„f (/Cnj.) BTII . Such a fc matrix can be established for a larger cluster, too. 

When expanded to lowest order in 6, the blocks of k may be diagonalized. It is then easily seen that 
they may be written as e~ s ( E< *~ flNa > with E a the eigenvalues of the Hamiltonian Eq. (Q3. In particular, in 
the two-particle sector, the eigenvalues 




,(.s) 



1 



V± 2 



^3e c + e a + I ± yj{e c - e a - I) 2 + 8V 2 j (25) 
correspond to the singlet states, and the three-fold degenerate eigenvalue 



= e c + e a + / 



(26) 
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corresponds to the triplet states. This exact evaluation of the partition function allows us to conclude that 
the functional integral representation of this SIAM, Eqs. (|9WTT1i. together with Eq. (f2Tb is correct. 
Should we use 



one encounters ill-defined integrals with the evaluation of V n ■ L an K an . Consequently this form of little 
use. Therefore, while various forms for y and z yield the correct fc-matrix, Eq. (T24t . they nevertheless need 
to be carefully verified. They are equivalent on the operator level, but their respective field representation 
involves a discrete time dependence and the time steps have to be tuned meticulously so as to render the 
functional integrals finite. 

Moreover, while there is no preferred form for y and z when computing the path integral exactly, 
their precise choice makes a difference on the level of the saddle-point evaluation ||29] [30l (for a recent 
comparison, see Refs. l2Tl[32l ). These seemingly formal properties ensure that the approach captures 
characteristic features of strongly correlated electrons as the suppression of the quasiparticle weight and 
the Mott-Hubbard/Brinkman-Rice transition ll33l to an insulating state at half filling with increasing on- 
site Coulomb interaction. In fact, the Kotliar-Ruckenstein approach has been impressively successful when 
compared to numerical simulations: ground state energies [29 J and charge structure factors show excellent 
agreement Il34l . 

3 Expectation values, correlation functions, and thermodynamics 

Having gauged away the phases of the slave boson fields leaves us with the amplitudes of the slave boson 
fields that are gauge-invariant quantities. In contrast to ordinary Bose fields, the expectation values of 
which typically vanish due to phase fluctuations, radial slave boson fields are generically characterized 
by finite expectation values that cannot be "nullified" by phase fluctuations. Indeed, they are intrinsically 
phase stiff. 

In order to illustrate the above let us determine (i?<r,i). According to the above it is given by a specific 
projection of the fermionic determinant, namely: 




(27) 



with e = + , instead of Eq. (|2TT >. leads us to the same conclusion. 
In contrast, when using 




(28) 




(29) 



with 
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fcf = ( ) ffi ( „ „ ] ( „ T I ffi ( ) © ( L c L a e 




-81 



and 



( 






° ) 




V o 






o ) 


ffi( 



(30) 



) © 



L c L a e- 51 ) © ( L 2 c L,e- 2SI ) © ( ) (31) 



Carrying out the algebra (at T = 0) yields: 

(i^l) = (l 2T/2 ; = W 1 

/ 4V 2 \ 

+ 1 = Uv 2 + <W 3 (32) 

\ (A - I) 2 + 8V 2 + (A - 7) V /(A - I) 2 + 8V 2 J ' ' 

where M is the number of electrons in the system, and A = e c — e a . Therefore (R a .i) arises as a specific 
projection of the fermionic determinant. It is non-vanishing for finite values of the particle number, V and 
A, thereby putting the above conjecture on firm grounds. 

Regarding the correlation functions they may be easily obtained from the above steps. For instance, the 
charge autocorrelation function reads: 



Z((R t:1 + R U i)(R t , m + R^ m )) = Tr((fct + h)k m ~ 2 {^ + k±)k N ~ m ) (33) 

Again, it is obtained as a specific projection of the fermionic determinant. This holds true for the spin 
autocorrelation function as well. 

With all the energy eigenvalues correctly reproduced, thermodynamics can be carried through. Re- 
markably, the specific heat CV=2 only depends on A — I. For example, for A — I < 0, it exhibits a 
Schottky-like shape for a two-level system, consisting of the levels and if the single particle state 
with e a is only thermally occupied and correlations play only a minor role. For A — I > the specific heat 
peak is controlled by (E^ — E^), it sharpens significantly, but also a weak high-temperature structure 

at T oc (E_l' — E^ ) is observed (upper right of Fig. [T). In the latter regime, correlations play a dominant 
role as the single particle state with e a is occupied. 

4 Summary 

Summarizing, we have revisited the Kotliar and Ruckenstein slave boson representation for the U = oo 
single impurity Anderson model. It backs on the introduction of three complex slave boson fields, subject 
to three constraints. The £7(1) x U(l) x £7(1) gauge symmetry of this representation allows to gauge 
away the phases of all three slave boson fields in the continuum limit. We then implemented a path in- 
tegral procedure involving radial slave boson fields defined for discrete time steps. The correctness of 
several Kotliar and Ruckenstein renormalization factors has been verified through the exact evaluation of 
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heat peak with A — I. 

the partition function and expectation values for a strongly interacting two-site model. In particular, the 
expectation value of the radial slave boson fields was shown to be finite, and does not refer to the introduc- 
tion of spurious Bose condensates. On the contrary they directly represent the density of empty or singly 
occupied sites. Correlation functions and thermodynamics can be correctly captured in this framework. 
Non-local Coulomb interaction can be taken into account as well, paving the way for calculations in the 
lattice case. For this purpose suitably modified ^-matrices, which enter the partition function Eq. (l22t . have 
to be identified. 
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